Measuring Probabilities which Violate "Reality" in a Bell Inequality Experiment 
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Bell inequalities limit measurable correlations between properties of two physical systems, under 
two assumptions: that these properties can be ascribed independent of measurement, and that influ- 
ences between them are bounded by light-speed propagation. Certain quantum entangled systems 
do not submit to this local realistic description, illustrated by copious experimental tests. Usually, to 
avoid quantum mechanical measurement back-action, different measurements are made on each of 
several subensembles of identical quantum systems. Here we show that weak measurements, which 
avoid back-action by weakly coupling the measurement device to the system, can be used to identi- 
cally measure all quantum systems in the ensemble, providing a direct empirical footing from which 
to reason about their properties. The quantum mechanical violation of Bell inequalities is explained 
by the negative weak-valued probabilities that we measure, compatible with an operational physical 
model which preserves locality at the expense of objective realism. 



Entanglement is the most striking property of quantum 
physics, providing a resource for quantum information 
technologies PUS] and leading to highly counterintuitive 
effects in the laboratory. Specifically, measurements on 
bipartite entangled systems demonstrate a contradiction 
with the notion that systems have properties that are lo- 
cally describable, independent of measurement, as first 
discussed by Einstein, Podolsky and Rosen (EPR) [6]. 
Following the proof that this notion — formalized as the 
joint assumption of locality and realism — leads to em- 
pirically testable Bell inequalities [7], many experiments 
(e.g. refs [SHE]) have demonstrated that at least one of 
these assumptions must be false, contrary to our strong 
classical intuition. 

Tests of Bell inequalities typically involve measuring 
two pairs of observables, one pair on each system. How- 
ever, because the observables in each pair are comple- 
mentary (e.g. spin projections ax — X and az = Z of 
a quantum bit), measurement back-action makes it im- 
possible to simultaneously measure all four combinations 
of observables on a given bipartite system. Rather, only 
one of the four combinations can be measured on each 
instance of the state. The inferences about the non-local- 
realism of the particles' descriptions come from the prob- 
abilities of outcomes on different sets of measurements on 
four different (albeit nominally equivalent) subensembles 
of quantum systems. 

Here we use weak measurements of polarization- 
entangled photons, and the formalism of weak values, to 
empirically determine joint probabilities for the simul- 
taneous values of complementary observables. In weak 
measurements, measurement back-action can be made 
arbitrarily small by reducing the strength of the coupling 
between the measurement probe and the system. Thus, it 
is possible to perform a weak measurement of one observ- 
able (Z, for example) followed by a strong measurement 
of a complementary observable (e.g. X) on the same par- 



ticle; the same technique can be utilized on the second 
particle (Figure [T]) . Although the results of weak mea- 
surements have a large variance due to the weak coupling, 
a good signal to noise ratio can be achieved by averaging 
over a sufficiently large ensemble. Because the same mea- 
surements can be applied to all members of the ensem- 
ble, this technique provides an observationally grounded 
method from which to infer properties of the quantum 
systems between preparation and the final strong mea- 
surements |15j . Our method is similar to that recently 
used in inferring the trajectories |16l I17| and wavefunc- 
tions [IS] of single photons in a double-slit interferometer. 
In the present case, our method allows us to examine the 
consistency of measurement results with a local but non- 
realistic description of the entangled state as introduced 
by Hofmann [IS], providing a new perspective on Bell 
inequalities. 

We gain new insights into quantum correlations and 
the foundations of quantum mechanics by applying 
weak measurements to the well-known Clauser-Horne- 
Shimony-Holt (CHSH) version [20 of Bell inequality. 
The CHSH inequality can be written as a bound on the 
expectation value of the CHSH parameter, which is a 
function of measurement result correlations between two 
systems. The CHSH parameter depends on the outcomes 
of two measurements on each of the systems: 

Schsh = (X + Z)P+(X-Z)Q, (1) 

where X, Z £ {±1} are the results of a pair of measure- 
ments on a system A (Alice's system), and P,Q £ {±1} 
similarly for a system B (Bob's system). Assuming local- 
ity, we can evaluate the expectation of this by summing 
over the 16 possible combinations of measurement out- 
comes, 

(Schsh) = K X + Z ^ P + ( X - Z ^ Pr t X ' Z - P ' 0]- 

X.Z,P,Q 

(2) 
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FIG. 1: Measurements in Bell inequality tests. A source pro- 
vides a pair of (possibly correlated) systems, one each to Alice 
and Bob, who each measure one out of two possible observ- 
ables. Traditionally (a) Alice and Bob each choose one or the 
other observable to measure strongly for each system. Using 
weak measurements (b) no such choice is necessary, as Al- 
ice and Bob may collect statistics of every observable with a 
fixed configuration. Doing so, Alice and Bob may infer prob- 
abilities for simultaneous outcomes that cannot be measured 
directly. Alternatively (c), Alice may perform weak measure- 
ments while Bob performs strong measurements. Configura- 
tions b and c generate identical results owing to the factor- 
ization of the CHSH parameter (see text). 



Notice that only one of (X + Z)P or (X - Z)Q con- 
tributes to each term of the sum because one must be 
zero. If the probability of each outcome is realistic, i.e. 
< Pr[X, Z,P,Q] < 1, it is simple to arrive at the CHSH 
inequality, |(5 C hsh)| < 2. 

The canonical example in which quantum mechanics 
violates this inequality is where Alice and Bob each pos- 
sess one part of the maximally entangled 2-qubit singlet 



state, I*" 



(|0) A ®|1> 



1 



>\0) B )/V2. The CHSH 



parameter becomes the operator 



'CHSH 



= (X + Z)®P+(X-Z)®Q. 



(3) 



To maximize violation of the CHSH inequality [21] , Alice 
chooses to measure from the Pauli bases X = |0) (1| + 
|1) (0|, Z = |0) (0| - |1) (1|, and Bob from the bases 
P= -(Z+A)/ v / 2andQ = (Z-X)/V2. A simple calcu- 
lation yields |(5 C hsh)| = |Tr[S CHS H |*~) (*~|]| = 2^2, 
violating the CHSH inequality by a factor of 

The key feature of weak measurements is that they 
minimize measurement back-action. This allows mea- 
surement results to be obtained for all four variables (X, 
Z, P, and Q) for the same two-qubit system simulta- 
neously. For example, weak measurements can be per- 
formed in the Z basis on system A and in the Q basis 
on system B, before strongly measuring the systems in 
the X and P bases respectively. Through this process 
one can determine |(Schsh)| using fixed measurements 
on every system in the ensemble. 

The measurement results have large variance due to 
the weak coupling of the measurement device to the sys- 
tem. However, the formalism of weak values allows us to 
extract a reliable, observationally justified value for the 



observables weakly measured on the ensemble. Within 
the overall ensemble of states with a given preparation, 
subensembles of systems can be defined by postselection, 
i.e. by the outcome of the final strong measurement. For 
each such subensemble, the average value of the results 
of the weak intermediate measurement, of some observ- 
able O in general, is known as the weak value |22j of 
O over that subensemble, and represents an empirically 
relevant way to talk about the value of O between prepa- 
ration and measurement. For qubits, every Hermitian 
operator is proportional to the identity plus a projector; 
the former has a weak value of one, and the latter has a 
weak value that can be termed a weak-valued probability, 
Pr w , for the eigenvalue associated with that projector. 
That is, for qubits, every measurement of a weak value is 
equivalent to measuring a weak- valued probability |23j . 

To see how this is applied to the CHSH inequality, first 
consider the weak measurement in the Z basis of a single 
qubit in a state \ip), followed by a strong measurement in 
X, which we can state as a postselection of the measure- 
ment eigenstate \<f>(x)) corresponding to the result x. For 
this postselected ensemble, the weak value, given by |16j 



Re 



(<Kx)\2\1>) 



(4) 



quantifies the expectation value of the results of weak 
measurements in the Z basis, as the measurement 
strength goes to zero. Weak values are unusual in that 
they can take values outside the spectrum of the measure- 
ment operator, but are nevertheless useful in increasing 
measurement precision where the resolution of the mea- 
suring device is otherwise the limiting factor 24-28] , and 
in resolving a number of quantum mechanical paradoxes 
(e.g. refs. [TBI ITT1 |2"5H3"5] ) and investigating macrorealism 
on one [3U [37] or two [35] systems. 

If we now consider the Z basis projectors, II 2 , the weak 
value j,(II z )^ can be interpreted as a weak- valued prob- 
ability of obtaining the outcome z, given the \tp) input 
state, conditional on finally postselecting the \<f>{x)) state 
indicating an outcome x in the X basis, i.e. Pr w [z|a;, ip] = 
0(II z )^. Because the back-action imparted onto the sys- 
tem by the weak measurement is negligible, it follows 
that the joint probability of obtaining both X — x and 
Z = z outcomes is therefore 



Pr w [x,z|^ = ,(H z );Pr[^)|^, 



(5) 



where Pr[^(x)|^] = \{(f)(x)\ip) | 2 . Weak measurements 
thereby allow us to ascertain (pseudo-) probabilities for 
outcomes that are not directly obtainable by typical mea- 
surements. 

Similarly, the joint probabilities for the two-qubit case 
can also be found: 

Pr w [x,z,p,q\il>]=t (ft, ® n 9 ); PtMx, p)\ V>] , (6) 
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FIG. 2: Experimental apparatus. Pairs of entangled photons 
are generated by a type-I bismuth borate (BiBO) 'sandwiched 
pair' source. One of the photons undergoes strong measure- 
ment, the other is coupled to single-mode optic fibre and di- 
rected to the weak measurement apparatus, implemented as 
a polarization interferometer. A tilted optical plate is placed 
in each interferometric arm, causing a polarization-dependent 
transverse shift of the spatial Gaussian mode of the photon, 
with polarization postselection following. The mode intensity 
is sampled by a scanning slit, and is collected into a multi- 
mode optical fibre (after passing through telescoping lenses, 
not shown). Processing the single-photon counting module 
(SPCM) detection rates completes the weak measurement. 



where H g represents projectors in the Q basis (applied 
here to Bob's qubit) corresponding to the outcome q, and 
the postselection state \4>{x,p)) depends on strong mea- 
surement outcomes of both Alice's (X = x) and Bob's 
(P = p) qubits. 

We can thus calculate the weak- valued joint probabil- 
ities of each of the 16 outcomes of measurements on the 
entangled state \^~)- They each take one of the four 
possible values: (2 + \/2)/16, \/2/16, (2 - \/2)/16, and 
— V2/16. Of particular note are the conditions for each 
of the two possible outcomes for the CHSH parameter, 
<Schsh = ±2. We find that the total probability for the 
positive outcome is (1 + v / 2)/2 ~ 1.207, while the prob- 
ability for the negative outcome is (1 — v / 2)/2 ~ —0.207. 
It is easy to see that the contribution to the expectation 
value, equation ([2]), is therefore positive in both cases, 
giving |(Schsh)| = 2\/2 as expected [T5] . 

These non- realistic (i.e. outside the range [0, 1]) prob- 
abilities may at first seem nonsensical. Indeed, they 
cannot be measured as relative frequencies in the 
laboratory — instead they must be inferred from weak 
measurement data, as we do below. They arise sim- 
ply from the fact that they stand for strong measure- 
ment results which cannot be physically obtained (events 
that cannot actually take place) because of measurement 
back-action. As shown previously pUJ [3S1 [33] j such non- 
realistic probabilities are nevertheless useful in probabil- 
ity accounting for physically realizable events. 

We now present an experimental demonstration of 
these non-realistic probabilities using photons. Firstly we 
note that if Alice could obtain both X and Z measure- 
ment results simultaneously, then by equation (nj) either 



P or Q of Bob's results would not matter for each shot. 
The only empirically relevant weak-valued probabilities 
are then either Pr w [cc, z,p] or Pr w [x, z, q]. We may there- 
fore simplify the experiment by performing weak mea- 
surement of only one photon of each pair (Alice's), rather 
than both simultaneously, with the probabilities for each 
outcome of the CHSH parameter remaining unchanged. 

Pairs of photons entangled in polarization and having 
high fidelity with | , F~) are generated via spontaneous 
parametric downconversion (using a type-I 'sandwiched 
pair' source [3D] : see the Appendix). One photon of the 
downconverted pair is assigned to Bob and immediately 
undergoes a strong measurement, either P or Q, by post- 
selection using a quarter- wave plate (QWP), half- wave 
plate (HWP), and polarizing beamsplitter (PBS). The 
other photon is assigned to Alice, coupled to a single- 
mode optical fibre, and guided to the weak measure- 
ment apparatus. Photons exit the single-mode fibre and 
pass through a collimating lens, resulting in a Gaussian 
mode. Photons then pass through a HWP at its op- 
tic axis. Weak measurements are achieved by engineer- 
ing a polarization-dependent displacement, Ar, signifi- 
cantly smaller than the width of the transverse Gaussian 
mode [H] . 

To achieve this displacement, we construct a partially 
spatially-mismatched polarization interferometer. A po- 
larizing beam displacer (PBD) separates the horizontal 
and vertical polarizations of the photon into two par- 
allel spatial modes. An optical plate is placed in each 
mode and tilted in approximately equal opposing direc- 
tions, causing small opposite displacements of the two 
modes following Snell's law. A HWP flips horizontal and 
vertical polarizations of the photon, the modes of which 
then (partially) recombine by the final PBD. Despite the 
interferometric modes not perfectly overlapping, a large 
Gaussian beam width relative to the size of the displace- 
ment ensures that decoherence (i.e. measurement back- 
action) is negligible. The photon is then postselected in 
the conjugate basis using a HWP at ±22.5° from its op- 
tic axis followed by a PBS. A slit, seated on a motorized 
translation stage, is scanned in the direction transverse 
to the beam, allowing us to sample the photon flux of 
the Gaussian profile at any point r within the range of 
the stage (25 mm). Integration of the photon flux (in 
coincidence with Bob's photons) allows us to determine 
the expectation of the weak polarization measurement. 
The joint probability of each polarization outcome of Al- 
ice's and Bob's photons can then be inferred (see the 
Appendix). Assigning horizontal and vertical polariza- 
tions to the Z = +1 and Z = — 1 states, respectively, 
we thus obtain weak measurements in the Z basis, with 
postselection in the X basis. 

The performance of the photon source is quantified 
by conducting quantum state tomography [12] with the 
slit removed and a QWP inserted immediately before 
the final HWP. With a two-photon state of measured 
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Bell-state fidelity 0.948 ± 0.002 and tangle g3] 0.841 ± 
0.004, we expect to be able to achieve a CHSH value of 
|(<Schsh)| ~ 2.67, sufficient to violate the CHSH inequal- 
ity. For this state, the experimentally determined weak- 
valued probabilities are 1.172±0.008 for Schsh = +2 and 
-0.171 ±0.002 for S CHS h = -2, resulting in | (Schsh) I = 
2.686 ±0.017 and violating the CHSH inequality by more 
than 40 standard deviations. 

Our choice of weakly measuring in the Z basis is 
arbitrary — we may instead measure X weakly before a 
strong Z measurement and obtain essentially identical 
expected joint probabilities. To demonstrate this, we in- 
troduce a change of basis of the measurements by set- 
ting the initial HWP in the weak measurement device 
to 22.5° from its optic axis. Performing the experi- 
ment in this condition we find a tangle of 0.857 ± 0.005, 
and probabilities of 1.147 ± 0.008 for S C hsh = + 2 and 
-0.140 ± 0.002 for S C hsh = -2. This gives a CHSH 
value of 2.574 ± 0.016. Despite the improved tangle for 
this state, sensitivity to imprecision in the manual set- 
ting of the change-of-basis HWP leads to slightly reduced 
magnitudes of the measured probabilities. 

With weak measurements in the Z basis, we also per- 
form the experiment for states of various tangle, obtained 
by tuning the source (see the Appendix). The results, 
shown in Figure [3^ a), demonstrate beyond- unit and neg- 
ative probabilities for states with high tangle (above 
3 - 2^2 ~ 0.17), with magnitudes decreasing with tan- 
gle. Violation of the CHSH bound is also seen for these 
states (Figure [3^b)). States with tangle below «0.17 do 
not exhibit these non-realistic probabilities. This corre- 
sponds to the CHSH value for these states passing below 
the bound of 2. (Any pure state with nonzero tangle can 
violate the CHSH inequality [43] . however this requires 
optimization of the measurement bases, which we do not 
do here.) 

For our apparatus, the results that we measure can 
be reproduced by a local hidden variable (LHV) the- 
ory; this is the case with any protocol in which there 
is no measurement choice on one side. However, the 
weak measurement formalism allows us to extract strong- 
measurement joint probabilities which cannot be repro- 
duced by any LHV theory. While standard CHSH tests 
may also be explained by recourse to negative joint prob- 
abilities (by assuming that equation ^ always holds and 
X, Z,P,Q £ {±1}), there is no a priori reason to expect 
such negative joint probabilities to have an experimental 
significance. In our weak measurement scheme, the neg- 
ative joint probabilities naturally arise as the outcomes 
of experiments probing the intermediate state of the sys- 
tem. 

Furthermore, the non-realistic probabilities arising 
from our weak measurement approach connect to a 
framework in which even entangled states can be de- 
scribed in a local way. Specifically, entangled states 
can be written as a statistical sum of local 'transient' 



e 0.5 




0.2 



0.4 0.6 
State tangle 



0.8 



2.5 



x 



1.5 



1 U 




0.2 



0.4 0.6 
State tangle 



FIG. 3: Experimental results for states of various tangle, a, 
Joint probabilities of the outcomes corresponding to Schsh = 
2 (upper curve) and Schsh = —2 (lower curve), b, The CHSH 
value | (Schsh) | calculated using the joint probabilities. The 
experimentally measured values (points) closely match the 
ideal theoretical values (lines) as the tangle of the state varies. 
Horizontal error bars represent the range of tangles measured 
via tomography taken before and after data collection. Ver- 
tical error bars span one standard deviation. 



states |15| . i.e. in the form pi — X)/ P r [/K]-^i/> where 
Rif is a separable matrix — a tensor product of matrices 
describing each qubit separately — as discussed by Hof- 
mann |15| 119) . The transient density matrices Ri / com- 
pletely represent the state as it transits from preparation 
in state pi to final measurement described by separa- 
ble two-qubit POVM elements {-?/}, yielding outcome / 
(reference |15jh They are given by 



Rif - 



F f pi + pjFf 



(7) 



Each satisfies two of the three conditions for an ordinary 
density matrix: it is Hermitian with unit trace. How- 
ever they may violate the third condition (positivity) 
by possessing negative eigenvalues. Nevertheless, these 
transient density operators provide empirically verifiable 
predictions of measurement results. 
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We used our strong measurement tomography of the 
input state p, to construct the various Ri / for our exper- 
iment, using equation Q. From these, it is possible to 
derive the joint measurement probabilities Pr w [ScHSH — 
+2] and Pr w [5cHSH — — 2], which we compared with the 
experimental joint probabilities found with weak mea- 
surements. From the Rif representation of the most 
entangled state in Figure [| (tangle 0.841), we found 
Pr w [5 CHS H = +2] = 1.168 ± 0.002 and Pr w [5cHSH = 
—2] = —0.168 ± 0.002, within two standard deviations 
of the weak measurement experimental results. We also 
found good agreement for all other tested states. 

As a result of equation |7|) , the weak values framework 
suggests an interpretation of Bell inequality violations in 
which the assumption of locality is retained and that of 
realism is discarded. This is because the measurement 
result / simply selects, out of the statistical distribution, 
one of the local subensembles Rif, which completely char- 
acterizes the coherences present in pi. Therefore, corre- 
lations in two-qubit measurement outcomes do not cor- 
respond to nonlocal action, but rather result from a clas- 
sical Bayesian updating of probabilities. On the other 
hand, the negative eigenvalues of the transient density 
matrix indicate that the system does not possess a real- 
istic state prior to the final measurement. 

By exploiting the power of weak measurements, we 
have presented an experimental violation of the CHSH 
inequality in terms of non-realistic weak-valued proba- 
bilities which align with a local interpretation of quan- 
tum mechanics. Our results point the way to further 
novel investigations of fundamental quantum phenomena 
in which typical strong measurements are replaced with 
weak measurements, possessing as little back-action as 
desired. As with the CHSH inequality shown here, this 
approach could offer deep insights into the most puzzling 
aspects of quantum mechanics. 

We thank Holger Hofmann for helpful discussions. 
This work was supported by the Australian Research 
Council. 
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Appendix 
Experimental joint probabilities 

Let th represent the centroid of the Gaussian pro- 
file associated with horizontally polarized photons ex- 
iting the interferometer, and ry represent the centroid 
of vertically polarized photons. Because ry > ra in 
our apparatus, it is convenient to write the displacement 
Ar = ry — Th ■ With horizontal polarization correspond- 
ing to Z = +1 and vertical polarization corresponding to 
Z = — 1, measuring the position of the Gaussian profile 
for a given photon state corresponds to a Z measurement. 
In typical spatial-mode qubit encodings, the modes are 
sufficiently far apart that the Guassian profiles are ap- 
proximately orthogonal, and the measurement outcomes 
thus correspond to strong measurements. The expecta- 
tion value for these measurements may be written 



(Z) = 



ru + r v -2r 
Ar 



p(r) dr, 



(A. 



where p(r) is the (Gaussian) probability density of de- 
tecting a photon at a position r. (For notational con- 
venience, we omit the explicit dependence on the initial 
state \ip).) As Ar is reduced towards zero, the measure- 
ment becomes weak and the uncertainty of individual 
outcomes increases due to the overlapping Gaussian dis- 
tributions for each outcome. In the absence of postselec- 
tion, however, the expectation value is the same in both 
the weak and strong measurement cases. Therefore, 



E 



zPr w [z] = r 
J —i 



Tji + ry — 2r 
Ar 



p{r) dr. (A.9) 



Using the fact that the two weak-valued probabilities in 



equation (A.9) must sum to one, it follows immediately 
that 



Pr w fZ 



1 

Ar 



(ry - r)p(r) dr, (A.10) 



and 



Pr w [Z = -l] = J- / (r~r H )p(r)dr. (A.ll) 

The postselection outcomes become extra conditions 
on equations (A. 10) and (A.ll). Supposing Bob mea- 



sures P, Alice X, then p(r) becomes p(r)Pr[x,p], and 
Pr w [z] becomes Pr w [z, x,p\. The joint probabilities can 
be calculated from experimental counting statistics by 



considering equations (A. 10 1 and ( A.ll ). For a finite slit 



width, the integral becomes a sum approximation over 
the range of measured positions as dr becomes Sr, 



Pr-[Z = +W]«i-£( 



(ry — r)p(r) Pr[x 1 p]5r, 

(A.12) 



and similarly for Pr w [Z = —l,x,p]. 

The value of p(r) Pr[x,p] cannot be measured directly, 
but must instead be estimated from photon detections. 
Let C(r,x,p) represent the count rate of these detec- 
tions. Then p(r) Pi[x,p]Sr = C(r,x,p)/CT, where Ct — 
J2 r x p C{r, x,p) is the total number of coincident photon 
detections over all the outcomes of the measurements in 
X and P. The joint probability can therefore be esti- 
mated by 



Pv w \Z 



-l,x,p] 



J2 r ( r v ~ r )C(r,x,p) 



(A.13) 



and similarly for Pr w [Z = —l,x,p], Pr w [Z = +l,x,q], 
and Pr w [Z = —l,x,q]. The probabilities of the positive 
and negative outcomes of the CHSH parameter are de- 
termined by taking the appropriate sums of these results. 



Source 

Photons pairs entangled in polarization are gener- 
ated using a type-I bismuth borate (BiBO) spontaneous 
parametric downconversion 'sandwiched pair' source [30] , 
pumped by 410 nm light from a mode-locked frequency- 
doubled Tksapphire laser. Two thin (0.6 mm) BiBO 
crystals, one cut for downconversion of horizontally po- 
larized incident light, the other cut for vertically po- 
larized incident light, are placed back-to-back. Pump- 
ing by diagonally polarized light induces coherent down- 
conversion from the two crystals resulting in the two- 
photon polarization state (\H) \H)+e lip \V) \V))/V%, for 
some constant phase ip. By changing the polarization 
of the pump light we can adjust the proportion sat- 
isfying the phase matching conditions of each crystal, 
thereby tuning the degree of entanglement generated. 
For a pump polarization angle 9, the resulting state is 
|^(0)) = wa.9\H) \H) +e^cose\V) \V) with a tangle of 
sin 2 26>. 

Each photon of the generated pair is assigned to either 
Alice or Bob. Before being coupled to a single-mode 
optical fibre, Alice's photon passes through a HWP set to 
its optic axis. This HWP is tilted around vertical such 
that it compensates for <p and achieves a state having 
high fidelity with \^~) at the measurement apparatus 
(in the condition of a diagonally polarized pump beam). 



Photon measurement 

The photon assigned to Alice undergoes weak mea- 
surement through polarization-dependent transverse dis- 
placement caused by tilted optical plates. For our appa- 
ratus it was logistically convenient to use quarter-wave 
plates set at their optic axes. Upon successful postselec- 
tion the photon then arrives at the scanning slit, after 



7 



which the photon passes through telescoping lenses (not 
shown in Figure [2]), a 3 nm full width half maximum in- 
terference filter, and is coupled into a multi-mode fibre 
guided to a single-photon counting module (SPCM). Sim- 
ilarly, following postselection, Bob's photon also passes 
through a 3 nm interference filter, and is coupled into 
a single-mode optical fibre guided to a SPCM. Alterna- 
tively to postselection, one could also perform strong two- 
outcome polarization measurements. Neither the polar- 
ization filtering method that we use, nor the slit-induced 
spatial filtering, are fundamental. 

Coincidence counting 

Detecting Bob's and Alice's photons within a coinci- 
dence window of ~3 ns helps ensure high-fidelity entan- 
gled two-photon states, however some accidental coin- 
cident detections remain. We estimate these accidental 
count rates by simultaneously recording detection events 
of the two SPCMs coincident when displaced in time by 



the pump pulsing period. All such detections thus arise 
from uncorrelated events. We subtract these accidental 
coincident counts from our results. 



Data collection 

We use a slit of width approximately 350 um, scanned 
over a 3.5 mm range in steps of 87.5 um, counting for 10 s 
at each step. This is performed for each of the eight total 
postselections for A, P, and Q, in turn. The process is 
repeated 70 times, ensuring that any drift in pump power 
is experienced equally (approximately) for each outcome. 
To calculate the joint probabilities it is necessary to char- 
acterize the centroid positions rjj and ry. This is done by 
additionally postselecting horizontal and vertical photon 
polarization states, making a total of ten measurement 
conditions in the experiment. We fit a Gaussian function 
to the count rates of each characterization — the centres 
of these fits define Vh and ry. 



